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The loop-current state discovered in the pseudogap phase of cuprates breaks time reversal sym- 
metry and lowers the point group symmetry of the crystal. The order parameter and the magnetic 
structure within each unit cell which is associated with it can be described by a toroidal moment 
parallel to the copper-oxide planes. We discuss lattice point group symmetry of the magnetic struc- 
ture. As an application, we discuss a few effects that necessarily accompany order parameter in 
the pseudogap phase. The magnitude estimated for these specific effects makes them hard to ob- 
serve because they rely on the small magnetic fields associated with the order parameter. Effects, 
associated with the electronic energies are much larger. Some of them have already been discussed. 

PACS numbers: 74.72.-h,75.25.+z,72.55.+s 



I. INTRODUCTION. 

Polarized neutron scattering experiments, dichroic an- 
gle resolved photo emission experiments and magnetiza- 
tion measurements^! 3 * 4 * 5 ! reveal the predicted line of 
phase transition^ in cuprates at pseudogap temperature 
T*(x). At the present time four classes of cuprates have 
been studied, all of which show consistency with the same 
form of order. The purpose of this paper is to present a 
symmetry analysis of the order parameter, which is neces- 
sary for further work. The neutron spin-flip intensity ob- 
served in the experiment in the pseudogap phase is con- 
sistent with the magnetic structure that transforms under 
the two-dimensional irreducible representation E u of the 
point group of the copper-oxide lattice, D^. The order 
parameter L that can be used to characterize the mag- 
netic structure of such symmetry is a polar time-reversal- 
odd vector parallel to the copper-oxide planes. Such a 
vector has been termed a "toroidal" or an "anapole" mo- 
ment in the literature. The relation of such an order pa- 
rameter to the magnetic structure and the loop-currents 
in the unit cell of cuprates is specified. 

Several experiments other than th ose which measure 
the order parameter directl ^ * 2 * 3 * 4 * 5 ! have already been 
proposed which lead to unusual effects due to the cou- 
pling of the order parameter to external probes. These in- 
clude several forms of dichroism in x-ray scattering^ and 
second harmonic generation in optical experiments^B. As 
an application of the symmetry analysis presented for the 
first time in this paper, we discuss few other physical phe- 
nomena that necessarily accompany the order parameter 
L. These are crystal lattice distortions which are second 
order in the order parameter, crystal distortion linear in 
an applied uniform field and linear in the order parame- 
ter, and magneto-electric phenomena in transport prop- 
erties. These are effects related to the small energies as- 
sociated with the magnetic field due to the order parame- 
ter. Much larger effects are associated with the electronic 
energies associated with the order parameter. This dis- 
tinction is similar to the difference between the exchange 



energies and energies of magnetic fields associated with 
spin-moment order. For example, the spin-splitting in 
ferromagnetic iron (ordered moment of « 2.5/is/unit- 
cell) is about 1 eV, while the magnetostrictive distor- 
tion due to the magnetic fields is a lattice distortion only 
about 1 part in 10 3 . 



II. LOOP ORDER PARAMETER: TOROIDAL 
MOMENT. 

The polarized neutron scattering experiments^! re- 
veal an elastic spin-flip intensity on top of the subset 
of Bragg peaks in pseudogap phase in cuprates. No new 
Bragg peaks appear; the magnetic structure responsible 
for spin-flip intensity docs not break translational invari- 
ance. However, it follows from the analysis of the neutron 
data that the symmetry of the observed magnetic struc- 
ture is lower than the point group symmetry of the lat- 
tice. The copper-oxide lattice has tetragonal symmetry 
4/mmm (D^h) (the effects of the small orthorhombic dis- 
tortion due to ordering in Cu-0 chains in some cuprates 
will be mentioned). In the pseudogap phase the symme- 
try is lowered to rrvmm or D 2 h(C2vr^- Such lowering 
of the symmetry follows if the magnetic structure trans- 
forms under irreducible representation E u of the point 
group of the lattice. For an order parameter in the pseu- 
dogap phase one can choose a polar in-plane (parallel 
to copper-oxide plane) vector L — (L y ,L x ) which is re- 
stricted to four crystalline directions because it trans- 
forms the same way, E u , under the operations of point 
group of the lattice. Order parameter L is time-reversal- 
odd because the magnetic structure which it represents 
is odd under time reversal. In the literature polar time- 
reversal-odd vector is called a toroidal moment as it is 
a symmetry of a magnetic field in a solenoid bent into a 
torua^ in the particle physics object of the same symme- 
try is known as anapole momentPA 

A natural way to relate the order parameter L to 
the pattern of the magnetization M(r) in the pseu- 
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dogap phase is as follows. Since the observed mag- 
netic structure retains lattice translation symmetries it 
is enough to consider magnetization M(r) within sin- 
gle unit cell. In general, the function M(r) can be de- 
composed into spatial harmonics of point group of the 
lattice (D 4h ); M(r) = Y.\ a c\ a M Xa (r) ■ Here har- 
monic M\ a (r) transforms under the irreducible repre- 
sentation A of the group D^h and a is internal index of 
the representation A if it is not one-dimensional; c\ a is 
a set of coefficients. The mathematical representation 
of the fact that the magnetic structure M(r) belongs 
to the irreducible representation E u is that it can be 
written as M(r) — c\M\(r) + C2M2(r) where Mi^(r) 
are two orthogonal harmonics that transform under two- 
dimensional irreducible representation E u . It can be seen 
that under all group operations a pair of vectors £ a =i,2 
defined by 



drM a (r) x r 



(1) 



unit cell 



also transforms under E u representation. The two vec- 
tors L a are parallel to the copper-oxide plane and are 
orthogonal to each other; we can use them as a basis to 
define a toroidal moment L = c\Li~\- C2L2 which is asso- 
ciated with the magnetic structure M(r) = c\M\[r) + 
C2M.2(v) in the unit cell. Conversely, if the magnetiza- 
tion M(r) does not have a component that transforms 
under E u , the integral in Eq. ([IJ vanishes. We also note 
that the out-of plane component of the integral in Eq. (JT|) , 
if non-zero, transforms under A2 U representation; it van- 
ishes if Afi 2(7") transforms under E u representation. We 
conclude that the vector L defined by 



drM(r) x r 



(2) 



unit cell 



can be used as an order parameter to characterize the 
magnetic structure in the pseudogap phase of cuprates. 

The simplest structure that has non-zero toroidal mo- 
ment is a pair of moments at finite offset; each moment is 
directed perpendicular to the line connecting them; the 
moments are of equal magnitude and opposite in direc- 
tion. However, function M(r) may have more elabo- 
rate structure; see the end of this section for discussion 
of the experimental situation. It is instructive to sep- 
arate M(r) into into planar and perpendicular compo- 
nents, M(r) = Mj_(r) + My(r); here M± = zM z and 
M|| = xM x + yM y and x,y is the basis in the plane 
parallel to the copper-oxide plane and z is in the direc- 
tion perpendicular to it. Fig. [2] schematically represents 
this decomposition with the red arrows understood as a 
magnetization direction M(r). Under the point group 
operations each component, Mi and My, transforms 
independently under its own E u representation. Conse- 
quently, each component may contribute independently 
to the toroidal moment L = (L x ,L y ) via Eq. p). We 
conclude that on symmetry grounds the magnetic struc- 



ture in the pseudogap phase can be an arbitrary combi- 
nation of , M± and M11 ; the precise balance between 
the two depends on the microscopic details and is decided 
experimentally, see discussion in the end of this section. 

Experimental data and theoretical calculations indi- 
cate that microscopic orbital currents within the unit 
cell are responsible for magnetic structure in pseudogap 
phase of cuprates. Here we discuss symmetry aspects of 
this relation. In general, magnetic structure in the crys- 
tal can be equivalently described in terms of the periodic 
magnetization function, M(r), or a pattern of periodic 
microscopic currents j m closed within the unit cell (loop 
currents). The two are related by 



j m = cVxM, 



(3) 



see Ref. [T3J In the absence of macroscopic currents, H = 
B + 4ttM = 0, the magnetization M(r) is proportional 
to the magnetic field within the unit cell, B — —4irM. 
The magnetization can be solved in terms of microscopic 
currents j m (r) by inverting equation j m = cVxM : 



M(r) 



1 

Aire 



dr 



i3mir') x (r-r') 



(4) 



where integral dr' is over the whole volume of the crystal. 
If the magnetization M(r) is a periodic function on the 
lattice, so is j m (r), and vice- versa. 

We can relate directly the toroidal order parameter 
L to the microscopic current distribution within unit 
cell. Current pattern is restricted to the unit cell, 
/unit ceil ^ r im( r ) = 0. Consider second moment of cur- 
rent distribution 



drr 3 m {r) . 



(5) 



unit cell 



Using definition of microscopic current in terms of mag- 
netization, jm — cVxM, we write this as 



/ drcr 2 VxM(r) . 

J unit cell 



(6) 



Using identity 

r 2 VxM = Vx(Mr 2 ) + 2M(r) x r (7) 

and the definition of the toroidal moment in Eq. ^ we 
find 

L = I dr(l/2c)r 2 j. m (r) - - [ ds x M(r)r 2 /2 . 

./unit cell v J 

(8) 

The second integral is over the surface of the unit cell, 
ds is surface clement. Using the fact that the surface is 
shared between adjacent unit-cells, the surface integral 
vanishes. Therefore, one can define a toroidal moment 

by 



L= I dr(l/2c)r 2 j m (r). 

i unit cell 



(9) 
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We conclude that the toroidal order parameter is pro- 
portional to the planar (parallel to copper-oxide plane) 
component of the second moment of microscopic currents 
within the unit cell. 

The magnitude of the spin-flip intensity observed in 
experiments for the most underdoped samples studied is 
consistent with a pair of magnetic moments ~ 0.1 fiB m 
the centroids of the two triangles; the centroids are lo- 
cated at distance xq ~ a/2 from the center of the unit cell 
(a is unit cell size in the CuO plane). The magnetic mo- 
ment M is given by a volume integral of magnetization, 
J drM(r). Integrating over half of unit cell (the full in- 
tegral over unit cell vanishes) we must obtain magnetic 
moment of order of 0.1/is; we conclude that the magne- 
tization in the unit cell is estimated as \M\ ~ 0.1/zs/a 3 
and the integral in Eq. Q is L ~ |M|a 4 ~ 0.1/iga. 




FIG. 1: (a) YBCO unit cell, black dots - copper atoms; open 
circles - oxygen atoms. Operations of the point group are 
indicated, (b) Representations of the point group can be con- 
structed on the space of patterns of currents on the inter- 
atomic links in the unit cell. Two diagrams indicate the pair 
of current patterns that transforms under inversion-odd two 
dimensional representation E u . The loop order parameter 
transforms under the same representation and can be repre- 
sented by these pictures. 

The simplest loop-current pattern that leads to mag- 
netic structure consistent with an in-plane toroidal mo- 
ment consists of a pair of planar current loops passing via 
planar oxygens and planar copper, see Fig. |TJb) . Such 
distribution of microscopic currents generates two oppo- 
site magnetic fluxes directed along z axis, see Fig. [2] 
as has been mentioned earlier, this is consistent with 
the pattern of the magnetization which has non-zero 
toroidal moment. However, polarized neutron scatter- 
ing experiments indicate a finite horizontal component 
of the magnetization M(r) in the unit cell^l. If the 
assumption is made in the analysis of the experiments 
that the structure factor for the in-plane and the out 
of plane components is the same, the magnitude of the 
horizontal and the vertical components are similar. To 
understand the presence of both horizontal and verti- 
cal components in the pattern of magnetization one has 
to consider also spin-orbit interactions and/or loop cur- 
rents involving the apical oxygen. In YBCO each copper- 
oxide plane is seeing non-centro-symmetric environment; 
consequently, symmetry allows a Dzialoshinskii-Moriya 
type interaction which in the presence of vertical com- 
ponent of the magnetization (planar loop currents) leads 
to a horizontal component of the magnetization due for 



example to an in-plane polarization of electron spin in 
the unit celP^. In mercury compound this mechanism is 
not allowed because each copper-oxide plane sees centro- 
symmetric environment. A horizontal component how- 
ever naturally app ears if one considers loop currents via 
apical oxygerpE2 as well as via planar copper. It is then 
reasonable that the major part of the horizontal com- 
ponent in YBCO is also due to loop-currents via apical 
oxygens. 




FIG. 2: Schematic representation of the magnetic flux distri- 
bution observed in the experiment. The diagram symbolically 
represents decomposition of the magnetic structure into hor- 
izontal and vertical components. The actual composition of 
the two into the physical order parameter depends on micro- 
scopic details and is not dictated by symmetry. 



III. LATTICE DISTORTIONS THAT 
ACCOMPANY LOOP-CURRENT ORDER. 

As mentioned earlier^, mmm symmetry does not al- 
low piezo-magnetism, i.e. no distortion changing sponta- 
neously the symmetry of the unit-cell to linear order in 
the order parameter is allowed. Here we discuss lattice 
deformations that accompany current-loop order to sec- 
ond order. To do so and for the other results derived in 
this paper, one must specify the irreducible representa- 
tions of the crystal symmetry in the absence of the loop- 
current order. The copper-oxide lattice has tetragonal 
symmetry 4/mmm (D^). The operations in the group 
fall into 10 equivalence classes: 7r/2-rotations around 
z axis C4, 7r-rotations around z axis C2, 7r-rotations 
around x axis C 2 , it- rotations around x + y axis C'2, 
7r/2-rotations around z axis followed by reflection S4, 
reflections <Jh,<Jv&d in the planes perpendicular to axes 
z,x,x + y respectively, spatial inversion i and identity 
operation E, see Fig. [TJa). Under symmetry operations 
in D^hi any physical object transforms under one of the 
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10 irreducible representations 



B 



E„ 



-i 2 2 

1 or z or x 
xy{x 2 - 



-y 
v 2 ) 

2 2 

ar — y z 
xy 

(zx, zy) 



Am 
A 2u 

B2u 
E,. 



xyz{x 2 — y 2 ) 

z 

xyz 
(x 2 -y 2 )z 



(10) 



The second column illustrates each irreducible represen- 
tations by a polynomial of the same symmetry Within 

the vector representation of the rotations in space is 
no more irreducible; instead, the the polar vector E and 
axial vector M break into irreducible representations of 

as follows : 



A 2g 
Eg 

A 2u 
E u 



M z 

(M y ,,M x ,) 
E z 

(Ey> , E X > ) 



(11) 



The basis x, y is chosen along the copper-oxide links. We 
use notation x',y' for a basis rotated by 45° with re- 
spect to x, y, i.e., x' = {x — y)/V% and y' = (x + y)/^/2. 
Representations E u and E g are two-dimensional repre- 
sentations, all other irreducible representations are one- 
dimensional. 

The term in the free energy that can couple loop or- 
der and lattice distortion has to be of at least a second 
order in the current-loop order parameter since L is time- 
reversal-odd; such a term has a structure LLu where u 
describes distortion. As a consequence, the lattice distor- 
tion is of a second order in loop-order parameter, u oc L 2 . 
To obtain the symmetry of allowed distortions we have 
to break bilinear products of L into irreducible represen- 
tations. Within D^h the product of two polar in-plane 
vectors (transforming under E u representation) breaks 
into four irreducible representations : 



Alg 


J 2 


+ Ll 


Alg 


Ly 


L x' ~ 


Big 


Ly 


L x ' 


Big 


T 2 





(12) 



Let us consider each of the terms above. One kind of dis- 
tortions that belongs to A\ g variety amounts to a change 
in the size of the cell and is obviously allowed to cou- 
ple to the square of the order parameter. Examples of 
the distortions in the other irreducible representations 
are shown in Fig [3] where black arrows represent atom 
displacements. The most interesting case is that corre- 
sponding to the third and fourth lines above: the free 
energy contains a term lB 2 g u B2g{L 2 i — L 2 ,) where ub 29 
is the distortion that belongs to B 2g irreducible repre- 
sentation and 7b 2 is a coupling constant; similarly for 
B\g. From the point of view of symmetry, this distortion 
is reminiscent of how finite polarization is generated in 
improper ferroelectric^^. 

In YBCO distortion of the A 2g symmetry is allowed 
in the presence of loop order and intrinsic orthorhom- 



bic distortion. At low enough temperatures the crys- 
tal structure of YBCO has intrinsic orthorhombic dis- 
tortion u(Big) of the Big symmetry due to the order- 
ing in the copper-oxigen chains (the lattice constants are 
a = 3.82A, b = 3.89A). Within D ih group the follow- 
ing relation holds Bi g B 2g = A 2g (this can be deduced 
from the character table). Therefore, there exists a term 
u(A 2g )u(Bi g )[L 2 , — L 2 ,] which will create a distortion 
u(A 2g ) of A 2g symmetry in the presence of finite or- 
thorombic lattice distortion u(Bi g ) and loop order pa- 
rameter L x i or L y '. Another effect of finite orthorhom- 
bic distortion u{Bi g ) in YBCO is to introduce a mixing 
between L y i and L x i via the term u(Bi g )L y i L x i . Along 
the same line of analysis one concludes that no distor- 
tion of A 2g symmetry is allowed in mercury compound 
in the pseudogap phase since its lattice has tetragonal 
symmetry. 

The magnitude of the distortions may be estimated 
from the free energy 



SF = ~f(u/a)L 2 + E{u/a) 2 /2 



(13) 



where u is distortion in a particular symmetry channel (a 
is in-plane lattice constant), L 2 stand for one of the loop- 
order bilinears in Eq. ( 12 1), and E is the elastic modulus 



(in energy units) for the distortion in the correspond- 
ing symmetry channel. Minimizing the Free energy we 
obtain u/a = jL 2 /E. It is quite difficult to get an esti- 
mate of the magnitude of the 7's directly. For an order 
of magnitude estimate, we proceed as follows: The Ai g 
distortion arises from the same symmetry arguments as 
the change in the volume of ferromagnetic substances due 
to their order and we may use experimental results ob- 
tained for them to get an idea of the order of magnitude 
of the effect to be expected. For example, detailed mea- 
surements are available for F^Q. The lattice constant 
changes from above the transition at over 800 degrees C 
to 300 degrees C due to ferromagnetism with a moment 
of about 2.5 /jb per atom by less than a part in 10 3 . With 
the ordered moment an order of magnitude smaller, we 
would therefore expect for similar coupling constant and 
Bulk modulus, a change in similar temperature region of 
only a part in 10 5 . This estimate is unlikely to be in- 
correct by more than an order of magnitude. One would 
expect that similar magnitude of distortion is to be ex- 
pected in the interesting case of the B 2g distortion and 
much smaller for the A 2g distortion since that must rely 
also on the orthorhombicity of the original structure. 



IV. DEFORMATIONS INDUCED BY 
MAGNETIC FIELD. 

We now consider possible unit cell distortions that can 
couple linearly to magnetic field H and linearly to the or- 
der parameter. Such distortion is allowed by symmetry, 
since the order parameter L is polar time-reversal-odd 
vector. So a term in the free-energy cx uLH is allowed 
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E„ 



FIG. 3: Inversion-even representations of the group D^h- 
Black arrows represent a polar vector, such as atom displace- 
ment; red arrow represent an axial vector, such as magneti- 
zation. 



with the lattice distortion u being odd under space in- 
version. To find the symmetry of all allowed distortions 
we must decompose the product of L and H into irre- 
ducible representations. For vector H = (H x >,H y ') in- 
plane (transforming under E g ), the product LH breaks 
into four irreducible representations 



Am 



Ly* Hy! ~~\~ Lr r ' H T. 



Ly> H T 



L x r H„, 



Ly' Hy! L X 'H X > 
Ly'H X ' + L X 'Hy' 



(14) 



For magnetic field out of plane H = H z (which trans- 
forms under A 2 g) the distortion u = (u x i , u y >) transforms 
under E u . Indeed, the product uL contains irreducible 
representation A 2g : 



A 



2fl 



Ly'lix' Lj_ 



V 



(15) 



We conclude that depending on the orientation of the 
magnetic field and the order parameter, distortion can 
appear in each of the five inversion-odd representations, 
see Fig|4j 

To estimate the magnitude of the distortion we con- 
sider free energy 






FIG. 4: Inversion-odd representations of the group Dah- 



introduces additional contribution to magnetic suscepti- 
bility via the term CL 2 M 2 which is obtained upon min- 
imization with respect to u/a; here C — f3 2 /E. The 
change in susceptibility in the presence of loop order 
is 5\ ~ ~X 2 CL 2 , see Ref [TBI for details. Both L and 
M can be measured in units of Bohr magnetons per 
unit cell, with L interpreted as staggered magnetiza- 
tion within unit cell; the convenient energy units are 
set with /ij^/A 3 ~ IK. In the experiment the suscep- 
tibility is measured 3mm 3 /Mol (or about 3 x 1CP 5 in 
dimensionless units for YBCO) at temperatures about 
pseudogap temperature. Since the loop order saturates 
rather quickly below the pseudogap temperature, we can 
attribute experimental value S\ ~ 0.1% to saturated loop 
order magnitude L s ~ 0.1 /its /Unit cell. We estimate 
C ~ ($X/X 2 )(1/L 2 ) ~ lO^-yUnit cell; the unit cell 
volume of YBCO is - 170l 3 . Using C = (3 2 /E with 
E ~ lOeU/Unit cell we estimate (3 ~ 5 x 10 6 in dimen- 
sionless units (where L 2 and M 2 arc in units of energy per 
volume of unit-cell) . We can now estimate the magnitude 
of the distortion. In the experiment the magnetization 
at external field of IT corresponds to magnetic moment 
M ~ 10 ~ 5 Hb /Unit cell. In this situation the distortion 
is estimated as u/a = (3LM/E ~ 1CT 6 . 



SF = (3(-)LM + E(-) 2 /2 
a a 

where u is distortion in one of the possible irreducible 
representations and LM is the corresponding bilinear. 
Minimization gives u/a — /3LM/E. We can estimate 
the magnitude of the coefficient (5 as well as typical dis- 
tortion magnitude using data from recent magnetization 
measurements-^. The cubic term of the form (3(-)LM 



Magnetic field dependence of elastic neutron 
scattering intensity. 



We now discuss a possibility of detecting lattice dis- 
tortions in the present of external magnetic field by the 
usual x-ray or neutron crystallographic techniques. The 
dominant contribution to elastic neutron scattering in- 
tensity comes from contact interaction between neutron 
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and atomic nucleus. One can analyze these starting with We rewrite these in a matrix form 



I N (q) cx \A N (q)f 

A n{0) = ^2 fa' 
R.a 



iq(R+r a ) 



(16) 



where q is (neutron) momentum transfer. r a is position 
of atom a within unit cell R; f a are nuclear scattering 
amplitudes. Let us assume that due to lattice distortion 
atoms shift to a new positions r a — > r a + u a within unit 
cell and ask whether there is a change in intensity to 
linear order in u a : 



5I N (q) = 6A^ N (q)A N (q) + c.c 



(17) 



For centro-symmetric crystal the sum in Eq. ( 16 1 evalu- 
ates to a real number (the nuclear amplitudes f a are real 
for q typical of slow neutrons), i.e., the amplitude A^{q) 
is real. To linear order in u a the change in the amplitude 
is 



5A N (q) = f a iqu a e 



iq(R+r a ) 



(18) 



R,a 



In loop-ordered phase, external magnetic field can pro- 
duce several types of distortions, all of which are odd 
under spatial inversion, see Eqs. (14 1 and (15 1. For the 
inversion-odd distortion u a the sum in the right-hand 



side of Eq. ( 18 1 is pure imaginary; intensity variation 
SIn(q) given by Eq. (17 1 vanishes. We conclude that 



the lattice distortion associated with external magnetic 
field, though non-zero, does not change intensity of the 
Bragg peaks in elastic neutron scattering to linear order 
in magnetic field. 



V. MAGNETO-ELECTRIC EFFECTS IN THE 
PSEUDOGAP PHASE 

Copper-oxide metal in the pseudogap phase exhibits 
magneto-electric effects. Magnetic field out of plane be- 
longs to A2 g irreducible representation. The product of 
an in-plane electric field E and loop order L has a com- 
ponent E x L that belongs to the same, A^ g irreducible 
representation (and similarly for out-of plane electric field 
which belongs to A 2u ) : 



A 2g 

A 2u 



E z 



E X L — -^y ' x ' -^x ' y ' 

H x L = Hy'L x i — H x /L y f 



(19) 



Corresponding to the two lines in this table, there are 
two terms in the free energy that couple electric and 
magnetic fields in the presence of loop order parameter 
L = (L x r,L y t) : 



SFme =Ai H z (E y 'L x > — E x 'L y 



A 2 E z (H y iL x i — H X lLyl) . 



(20) 



[H x . 











—\2Lyi 




E x i 


Hyl H z ] 








\2L X ' 




Ey, 




— \\Lyl 


XiL x / 







E z 



(21) 



The two coupling constants Ai^ have different magnitude 
since they describe coupling in two different irreducible 
representations. Assuming a single domain phase, say 
L = (L x i,0), we obtain a magneto-electric tensor in the 
symmetry-broken phase of the form : 




a 2 
ai 



(22) 



where ax,2 — X\.2L X . General symmetry consideration 
allow axial time-reversal-odd tensor of the second rank 
of this form in the crystal class mmm, see Ref. 1151 

We now discuss the magnitude of the magneto-electric 
susceptibility a in cuprates that is associated with the 
magnetic structure in the pseudogap phase. In spin sys- 
tems the magneto-electric coefficient is limited by the 
magnitude of the spin-orbit coupling^ a typical value 
a ~ 10~ 4 to 10~ 3 (dimcnsionlcss in cgs units) is several 
orders of magnitude smaller than the upper bound set by 
thermodynamics^!. In the cuprates the magnetic struc- 
ture in the pseudogap phase has an orbital origin. How- 
ever due to finite resistivity of cuprates in the pseudogap 
phase the observation of the magneto-electric phenom- 
ena may not be possible. To estimate the magnitude of 
the magneto-electric coefficient a we start with thermo- 
dynamic identity (see, e.g., Ref. [TU)l 



d 2 F 



dMi 



dR. 



dE.dHj dEi dHj 



(23) 



We use the last equality, i.e., we apply magnetic field H 
and estimate the induced polarization. It has been esti- 
mated earlier that in the experiments in Ref. QO the dis- 
placement of copper and oxygen atoms relative to each 
other is of order u/a ~ 10~ 6 at external field of IT; the 
symmetry of the displacement has been discussed earlier 
in the paper. The polarization which corresponds to such 
displacement is estimated as P ~ (e/a 2 )(u/a) <~ 10~ 3 esu 
(here a is the ab plane unit cell size ~ 4A; we have as- 
sumed that copper and oxygen in the copper-oxide plane 
have charge ±e). We find that magneto-electric sus- 
ceptibility is a — P/H ~ 10~ 5 , (which is dimension- 
less in cgs units). This corresponds to about 3 x 10~ 8 
Gauss/ (Volt /cm). This value could be observed exper- 
imentally in an insulator but not in the metallic pseu- 
dogap state of the cuprates where final voltage arises in 
the bulk only due to finite resistivity. We should men- 
tion for completeness that the orbital ordered phase also 
supports a zero-field Hall effect. But for a field applied 
of 10 volts/cm, the Hall effect would correspond only to 
that due to a magnetic field of about 10~ 6 Gauss. 
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Not all effects accompanying the observed orbital order 
are so small. For example, x-ray dichroism is predicted 
to be observable^. So is the predicted second harmonic 
generation^. All the effects investigated here come due 
to the periodic magnetic fields generated by the orbital- 
ordered phase, which are always small. The effective 
electronic energy of such phases is two to three orders 
of magnitudes larger than the magnetic energies. This is 
similar to the difference between exchange energies and 
magnetic field energies in magnetic order due to spin- 
moments. 
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